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IMPORTANT ANNOUNCEMENT 


The increase in the output of papers in the field of mathematical 
biophysics makes it difficult to insure prompt publication without an 
increase in the size of the journal. Therefore, the Bulletin of Mathe- 
matical Biophysics inaugurates the following service: 

Upon acceptance of a paper, the Editor, if necessary, will ask the 
author to shorten the paper to an extent dictated by the requirements 
of a reasonably prompt publication. The shortening should in no 
case reduce the paper to a mere abstract. Such a shortened paper will 
be published within six months or less. 

The unabbreviated original manuscript will be kept on file at 
the editorial office. Any person desiring to avail himself of the com- 
plete manuscript, may obtain promptly a microfilm copy of the latter, 
at the cost of 1¢ per page plus postage, by applying to the Editorial 
Office, 5822 Drexel Avenue, Chicago, Illinois. 

All papers in the Bulletin which have been thus shortened, will 
be marked at the end by the symbol MF, followed by a figure, indicat- 
ing the number of doublespaced typewritten pages of the unabbrevi- 
ated manuscript. 
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A MATHEMATICAL ANALYSIS OF ELONGATION 
AND CONSTRICTION IN CELL DIVISION 


H. D. LANDAHL 
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An equation for the rate of elongation of a dividing egg is integrated 
and generalized. The rates of elongation and constriction of a number of 
eggs under various conditions are analyzed and compared with the theo- 
retical predictions. The theory accounts rather well for a large body of 
data on elongation and constriction. The general shapes of the elonga- 
tion and constriction curves are predicted and the orders of magnitude of 
the parameters are satisfactory. One of the parameters for the elonga- 
tion curves is related theoretically to the parameter of the constriction 
curves, and the correct order of magnitude is obtained if one parameter 
is predicted from the other. 


The factors involved in cytokinesis are many and varied and the 
influence of various chemicals and physical changes on the process 
have been studied extensively. Several theories have been advanced 
as to the source of the mechanical forces involved in the process of 
deforming the cell (cf. Wilson, 1925). Included in these hypotheses 
are such factors as changes in form of the astral figure (Chambers, 
1924; Gray, 1931), production of a coagulating principle (Burrows, 
1927), variation in surface forces over the cell surface (Spek, 1918), 
forces of an electrical nature, changes in permeability (Lillie, 1918), 
changes in viscosity due to sol-gel transformations (Chambers, 1938; 
Marsland, 1939), equatorial intrusion and polar pseudopodial move- 
ments (Schechtman, 1937), osmotic forces and drag forces due to 
diffusion of metabolites (Rashevsky, 1938). 

It is possible to divide oil drops by reducing the oil—water ten- 
sion on two sides with a resulting elongation in the direction of mini- 
mum tension (Spek, 1918; Wilson, 1925). The possible effect of sur- 
face forces varying axially has been considered theoretically by G. 
Young (1939a) for the case of a true surface tension with a variation 
of the tension along the surface. Depending on how the tangential 
stress is sustained, one can obtain either shortening or elongation for 
a given case. The effect of surface forces will be discussed in greater 
detail subsequently. Calculations by N. Rashevsky (1938) and R. 
Williamson (1939) indicate that the effect of electrical forces would 
be small compared to those of diffusion. In the above mentioned hy- 


45 


46 MATHEMATICAL ANALYSIS OF CELL DIVISION 


potheses, except the last, little attempt has been made to predict any 
quantitative changes or to develop any mathematical theory of the 
effect of the various factors. Just how the forces arise is not clear in 
most cases. This in no way invalidates them, but makes difficult a 
mathematical treatment of the theories, at least until some investiga- 
tion is made concerning the manner in which the forces arise. How- 
ever, a considerable number of facts can be brought together with 
their aid. Theoretical investigation along these lines is definitely sug- 
gested. In the present paper we shall limit our discussion to the case 
in which only diffusion forces, osmotic forces and surface tension act, 
and in which viscosity opposes motion. ; 

The nature of diffusion drag forces has been discussed in a num- 
ber of papers (Rashevsky, 1938; Young, 1938; Landahl, 1942). As 
part of a program to investigate theoretically possible factors in- 
volved in cell division, this type of force has been considered (Rashev- 
sky, 1938; Young, 1939), and certain predictions regarding elongation 
and constriction have been obtained. It is our purpose to extend some 
of these results, compare the predictions with experimental data, and 
to discuss the implications of the restrictions thus imposed on the 
parameters. 


A COMPARISON OF THE ELONGATION EQUATION 
WITH EXPERIMENTAL DATA 


An equation for the elongation of freely dividing cells as a func- 
tion of time. We start with the expression for the rate of elongation 
given by N. Rashevsky (1940, chap. iii, eq. 44): 
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The parameters 27, and 2r, are the length and width of a cell 
whose diameter is 27, in the spherical shape; y is the surface tension, 
the viscosity of the interior of the cell, t the time, D, and D, the ef- 
fective external and internal diffusion coefficients, q the effective rate 
of production of metabolites of average molecular weight M, uw the 
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fractional volume of non-solvent particles on which the diffusion 
forces act, T the temperature and R the gas constant (Rashevsky, 
1940). In view of subsequent work (Landahl, 1942), the expression 
for K is actually more complex than given by (2), but about of the 
same order of magnitude. 


Equation (1) is not expected to hold for too great elongations. 
Thus let 


S(T — Fy) / To. (4) 
The volume is practically constant, and is given by 


ieee : nts — : wy. (5) 


From (4) and (5) we have 7; and r, in terms of ¢. Introducing these 
values into (1), expanding and retaining only the first three terms of 
the expansions, and rearranging, we obtain 
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Starting with a more complete equation (Rashevsky, 1940, chap. 
iii, eq. 14) including the permeability, h, so that osmotic forces are 
taken into account, we obtain as above an expression of the form (10). 
The parameters corresponding to A,B and C are now rather complex, 


and are of the form 
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(14) 


and the factors G, and G, depend upon the value of b as defined above, 
as well as upon the extent to which 37, exceeds D.(2 — 3u). The ac- 
tual expressions for G, and G, are very complex and so are not given 
here. 

As uw would not exceed about one-third, and is probably not less 
than one-twentieth, the requirements that the K’ > 0 is that 7h be 
larger than D, by a factor from } to 10 depending on uw. If 7h > D., it 
follows from the expression for A, the total diffusion resistance of a 
spherical cell (Rashevsky, 1940, chap. i, eq. 21 and 28), that the con- 
tribution to the resistance due to the external medium is greater than 
that due to the plasma membrane. This is not so for relatively imper- 
meable substances and in particular is not so for oxygen if one uses 
the known values of D., 7) and the estimated value for h (Landahl, 
1937; Rashevsky, 1940, chap. ii) for eggs with outer membranes. 
However, the permeability without the outer membrane may be con- 
siderably larger. If it is larger by .a factor of about ten, then we would 
have 7,h > D.. For many metabolites not enough information is 
known and thus we cannot arrive at any conclusions. The restriction 
is an important one, and must be considered in greater detail even- 
tually, but we shall not do so here. If mh < D., but if gq < 0, division 
is still possible for certain cell sizes. This case has been discussed by 
N. Rashevsky (1940, chap. iii). 

The effect of the variation of production rate with concentration 
on the elongation equation. In the development of the elongation equa- 
tion, it has been assumed that q was essentially a constant. We may 
inquire as to the effect of a g which depends upon the concentration. 
One might wish to consider the case in which the production, q, is 
proportional to the average concentration, or the case in which there 
is a constant production but a resynthesis going on proportional to 
the concentration. These cases can be summarized in the expression 
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d=Q+ac. (15) 


For the first mentioned case g; = 0 , while for the second case a < 0. 

Using the approximation method (Rashevsky, 1940), we may de- 
rive the expression for elongation. The expression is cumbersome, 
but if we proceed as above, we may obtain an expression for the rate 
as a parabola (see next section) in which instead of A we use the ex- 
pression A’ given by 


OF CLD) 
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A and B are again given by (7) and (8) but instead of q in K of equa- 
tion (2) we write q. + ac,, where c is the external concentration. 
The third term in the brackets of (17) is less than the second, since 
D. > D; and B >A. Thus, A’ and B’ differ from A and B by the fac- 
tor in the brackets of (16). Asq=q +ac,.,a< q/c fora>0O. For 
Arbacwa eggs the factor may be written approximately as one plus a 
quantity of the order of, or less than, qg/(be,.). Unless 6 and ¢ are 
taken to be smaller than would ordinarily be expected, this latter term 
is smaller than one, and A’ and B’ are approximately the same as A 
and B respectively. On the other hand, if the value of q/(bc,) is larg- 
er than one, A and B would be changed by about the same amount. If 
this term is considerably greater than one, the approximation method 
breaks down seriously. One could perhaps allow for a factor of the 
order of ten. 

If a is negative, A’ and B’ may be reduced in magnitude. If the 
factor q/(be.) is small, we may disregard it and we may drop the 
primes. If the factor approaches unity, then A’ may be appreciably 
less than A , and B’ might become negative. This is a rather special 
case. If the factor exceeds unity, A’ is negative. But we are only in- 
terested in A’ greater than zero, since this is the case for dividing 
eggs. 
From the above results, we see that the difference between a con- 
stant rate of production, q, and a q proportional to concentration or 
a q which increases or decreases linearly with concentration, is slight 
as far as its effect on the slope of the elongation curve is concerned. 
The parameters may be changed by a fairly large factor, but the ratio 
of A to B will probably not be affected appreciably. 
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The effect on the elongation curve of a force proportional to the 
gradient of the concentration compared with that due to the gradient 
of the logarithm of the concentration. In a previous paper (1942) it 
was seen that the force might be proportional to the gradient of the 
logarithm of the concentration c = mN, rather than to gradient c. 
The result is that, if the force is derived from a potential, the latter 
is proportional to log ¢. It is the difference in potentials at the “ends” 
and “sides”, log c, — log ¢., that enters into the elongation equation 
(Rashevsky, 1940, p. 43, for h = o). But this is log (¢:/c.). Only 
small differences of concentration are required to produce a consider- 
able effect. Thus, c, is of the order of c, and hence log (¢:/¢2), which 
is equal to log [1 + (¢; — ¢2)/c2.], is approximately (c, — C2)/t. 
Since the difference c, — ¢, varies much more than ¢,, we may treat 
€, as a constant for a first approximation. We then have the same re- 
sult that would be obtained if the force were proportional to gradient 
ce. Thus whether the force is proportional to gradient ¢ or to gradient 
log ¢, the resulting expression for elongation is the same to a first 
approximation. 

A general comparison of the equation of elongation with experi- 
mental data. Timed photographs of dividing Arbacia eggs during di- 
vison were made by R. Buchsbaum and R. Williamson (1942) for the 
purpose of testing some of the theoretical conclusions which follow 
from the diffusion drag force hypothesis. Some measurements of 
elongation of Arbacia eggs are given by L. Churney (1936, 1940), but 
no constriction measurements were made simultaneously as would be 
necessary for our discussion. We shall confine our considerations to 
the data first mentioned. In most cases the eggs were denuded of the 
fertilization membranes. These cells without membranes elongate 
markedly during the process, taking on a dumb-bell shape. It is this 
process of elongation which we are considering. After elongation 
ceases, the egg slowly shortens to some extent. For the present this 
effect will not be considered. 

The rate of change of the relative elongation with time is given 
by equation (6). The parameter A must be positive. A negative A 
would correspond to a case in which a cell became shortened before 
beginning to elongate. Since this is not observed, A must be positive. 
Then for A > 0, it follows that C > 0. The parameter B is positive 
if b = 1/8, and negative if b = 1. It changes sign at some point in 
this interval depending on the values of A and K. If B = 0 and 
C > B, we may ignore the last term and the expression is a parabola. 
On the other hand if B should be very much smaller than C we may 
drop that term and we have the rate given by the linear and cubic 
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terms. If B and C must both be taken into account, one will generally 
have a curve whose shape is intermediate. If one normalizes the ex- 
pression so that the maximum value of the rate, é'm , is unity, and the 
value of ¢, em, at which e’ returns to zero is also unity, then the 
Shapes of the curves may be readily compared. Designating the nor- 
malized variables as « and «*, the expressions corresponding to B = 0 
and C = 0 are respectively 


ets ook | 
e== SR Ou (18) 
2 
e* = 4e*(1 — «*), (19) 


These expressions have no arbitrary parameters and enable us 
therefore to compare data obtained on different cells. Thus, if the 
curves for these expressions are graphed, one should expect that any 
experimental data, to which they are to be compared, should lie be- 
tween them if there were no error of measurement. In Figure 1, these 
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FIGURE 1 


expressions are graphed, the former expression being shown as a 
broken line, the latter as a solid curve. To test this expectation, the 
measurements of the slopes of the elongation curves were made from 
each of thirty-four eggs under various conditions. These values were 
normalized in the sense used above and plotted in Figure 1. A discus- 
sion of the specific cases will be given later. The various symbols used 
to represent the measurements are explained in the figure. For each 
curve, the slopes were determined at about six to ten points. The scat- 
ter in the measurements is rather large as may be expected from meas- 
urements of slopes. However, within the range of error, it is clear that 
the average trend of the data lies approximately within the range of 
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the theoretical curves. The absence of arbitrary parameters in the 
theoretical curves increases the significance of the agreement. To the 
extent of the agreement we may say that the form of the experimental 
curve is predicted correctly by the theory. he 

It will be noted that the data scatter about the origin. The origin 
is determined in each case by the initial radius 7) as determined by 
proper average of several measurements of cell diameters prior to 
appreciable elongation. Some cells appear slightly elongated for 
some time prior to beginning of elongation for division while others 
are slightly compressed. A few showed a slight shortening before be- 
ginning to elongate. This may be due in part to changing in position 
of the cells as discussed by L. Churney (1936). If this is the case, 
the variation caused should be randomly distributed, and one would 
still be justified in using the initial trend. This effect is only appre- 
ciable near the origin. It should also be pointed out that the approxi- 
mations made in deriving the original expression (1) begin to break 
down at about the maximum of the curve. Without a more detailed 
investigation, one cannot expect agreement in the latter part of the 
curve, but it may be considered as an extrapolation. 

The fact that a parabola is a fairly good approximation to the 
experimental curve as shown in Figure 1 suggests that one might 
use only the first two terms of (6). This could be done only if the 
values of the parameters are not significantly altered. If one first as- - 
sumes C to be zero, and determines A and B from a comparison with 
the experimental curve, the parameters K and y/y are obtained from 
equation (7) and (8). Similarly one can obtain these parameters by 
assuming B = 0. It is found that essentially the same values are 
found for K and y/n in each case. The value of K obtained in the 
second case varies only by a factor 4 when b varies from 0 to 2, while 
y/n varies less than about 1/2 , unless b becomes large enough to make 
B negative when this can occur. But then C must necessarily play the 
dominant role and we thus obtain essentially the same values for K 
and y»/y, unless B takes on so large a negative value that a much 
larger value is obtained for C from the complete expression than 
would be the case when B is set equal to zero. We shall thus for sim- 
plicity drop the Ce? term as well as b in most of the further discus- 
sions. Instead of equation (10) we have 


e=A/(B + AeAltt)) , (20) 


where ¢, is an integration constant. A and B are given by expressions 
(7) and (8) so that we have, solving for K and v/n , 


K=7r,(A + 4B)/2, (21) 
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y/n = 2r,.(4B — A)/3. (22) 


E'ggs with membranes. In the process of division, eggs with mem- 
branes elongate about twelve per cent (cf. Churney, 1936). The mem- 
brane appears to act only as a restraining sac which is elastic, and 
can be considered as outside the cell proper. One should then prob- 
ably replace the surface tension y by a more complex expression. We 
yom however, evaluate the parameters on the basis of the same equa- 

ion. 

The curves of elongation of five eggs at first cleavage with ferti- 
lization membranes intact were plotted and the parameters A and B 
determined. One of these eggs had been treated in 10-*M KCN for 
twenty minutes and had been shaken. A second egg had also been 
shaken. The relative error is large because the amount of elongation 
is small, but the ranges of values for A and B were not large. The 
parameters A and B as well as the time from fertilization to first di- 
vision from the KCN treated egg did not take on values outside of 
the range of the other eggs. The average value of A for the other eggs 
was 0.025 + 5 sec. while the value of B was 0.22 = 8 sec. -?. No val- 
ues varied by more than fifty per cent of these averages. From these 
values, we find K = 1.7 X 10° cm. sec. and y/7 = 2.1 X 10° cm. sec.7. 
Further discussion of these values will be made subsequently. The 
average value of maximum elongation, ¢«,,, was é, — .12 + 1, while 
the average of the maximum slope, (de/dt) m = é'm, Was &m = 3.7 = 
9 X 10~ sec?. The average value of the time from fertilization to 
beginning of division t’,, as measured by the intercept of the line of 
maximum slope with the line «= 0, was found to be t’) = 48 = 3 min. 
The value t’, is more easily determined than t, of equation (20), and 
is about two minutes less than the latter, the latter being approxi- 
mately the time to the inflection point. The average initial diameter 
was 2r, = 74 + 2 X 10+ cm. The number after the + sign in each 
case gives the mean variation about the average. In this case it was 
estimated on the basis of values from the five eggs. 

Eggs with membranes removed by shaking. Sixteen elongation 
curves of eggs with the fertilization membranes removed by shaking 
were selected as being suitable to compare with the theoretical equa- 
tion (20). The parameters were thus found to have the values 
A = 0.024 + 4 sec. and B = 0.059 + 9 sec... From these values we 
obtain K = 4.5 X 10+ em. sec. and y/y = 4.9 X 10-* cm. sec.*. Other 
average values are: &m — .40 = 5, &'m = 2.3 Eo) 10-8 sect hg = 
46 + 4 min., 27, = 69.5 + 15 X 10% cm. Comparing these values with 
those from eggs with membranes, we see that not only the extent of 
elongation is greater without membranes, but the maximum rate is 
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about six times as large. The values t’, and 7 are essentially the 
same, as is also the parameter A. But B is much smaller, so that K 
and y/7 are both smaller, especially y/n. A sample curve is shown in 
Figure 2a. The curve of constriction will be discussed later. In this 
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case there are more points than there are for most curves. The para- 
meters are A — .029 sec.-?, B = .069 sec."?, t = 46.1 min., and 27) = 
75.5 X 10-* cm. This value is obtained by extrapolating the slope 
curve to zero and is two per cent higher than that obtained by averag- 
ing several diameters. 

Eggs with membranes removed by KCl treatment. Two elonga- 
tion curves from eggs with membranes removed by treatment with 
KCl were used to obtain the various parameters. The eggs did not 
show any hyaline material on the surface. The values of the para- 
meters are: A = .019 sec.*, B = .082 sec.-1, &, = .58, &'m = 2.7 X 10° 
sec.-4; fp = 44.7 min., 275 =.73.0:X 10% em: and A —.021 sec; Bo 
3029 ‘see.*): en == 71, eg 4 OL seceh t = 40.8 min, 2% oe 
X 10-* cm., giving the average values K = 2.6 X 10-* cm. sec.-? and 
y/n = 2.4 X 10+ em. sec.*. The latter curve is shown in Figure 2b, 
where ¢, = 42.9 min. Comparing with the values from eggs with mem- 
branes removed by shaking, we find that the extent of elongation is 
significantly increased, that the maximum rate probably increased, 
that A remained essentially unchanged, and B decreased about one- 
half, so that K and y/n decreased, especially the latter. 

Eggs with membranes removed by shaking and treated with 
KCN. Twelve elongation curves from eggs without membranes and 
treated with KCN were used to determine various parameters. There 
is some variation in both the concentration of KCN and time of treat- 
ment, but we shall treat the group as a whole. In most cases, the 
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concentration was between 10+ and 10-* M and the eggs remained in 
the solution from a short time after insemination. The average val- 
ues of the various parameters obtained are: A — .020 + 3 sec.-}, 
B= .048 + 9 sec, e, = .44 + 2, &'m = 2.2 + 4 X 10° sec, tp = 57 
= 7 min., 279 = 68.4 + 7 X 10+ cm., K = 3.6 X 10 cm. sec, and 
y/y = 3.9 X 10“ em. sec. Comparing these values with those from 
eggs treated the same but without KCN we find that the effect of the 
KCN is noticeable in the increase in t’,. The other differences are not 
at all significant. 

Eggs with membranes removed by shaking, second cleavage. 
Three cells with membranes removed by shaking were followed into 
second cleavage. From six curves, the following values were found: 
Agel +13) Secit.B == )071= 12 sec, a, = 85) 4 One een 8 
+ 3 X 10° sec.+, t’) = 72 + 2 min. from insemination, 27, = 58 = 1 
xX 10* cm., K = 4.5 X 10 em. sec.+, and y/y7 = 5 X 10+ em. sec.-. 
The values of 27, could not be measured with much accuracy, because 
the cells were not sufficiently regular. The values used here were de- 
termined by extrapolating the curves, and are perhaps ten per cent 
too high. Comparing the values of the parameters of second cleavage 
with those of first cleavage, we find the values are almost identi- 
cal, or not significantly different except for t’, and 27, which are neces- 
sarily different. The value for 27, for the half cells should be 0.8 that 
of the whole cell. A sample curve is shown in Figure 2c, where A = 
.025 sec.1, B = .065 sec.+, = 72.2 min., and 27, = 58 X 10+ cm. 

Discussion. In the early stages of the development of a theory, 
one does not expect too much agreement with experiment. When pre- 
dictions can be made, they should be tested. If the prediction is veri- 
fied by experiment, further development of the theory is suggested. 
If the prediction is not borne out, it is necessary to interpret the 
meaning of the failure on the basis of the theory, modify the theory, or 
suggest how other factors, which have not been considered may be 
responsible. It is essential that the order of magnitude of the para- 
meters be in agreement with values obtained by substituting plausible 
values of the more basic constants of which they are made up. The 
expression of K is given in equation (1). As discussed previously 
(Landahl, 1942), the expression for K is more complex, but the order 
of magnitude is probably not much changed. If we take u ~ 10°, 
M ~ 102, D. ~ 5 X 10%, and 7 ~ 107, then since K ~ 5 X 10%, we 
find that it is sufficient that q be of the order of 10° gm. cm.* sec.”. 
The parameter q is to be considered an effective rate of production, 
and thus may be of the order of, or less than the q values of the most 
rapid processes, those of respiration. This value is well below those 
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for respiration (Rashevsky, 1940, chap. ii). The value of K is thus 
found to be satisfactory. A possible increase in A and B, discussed in 
the preceding section, would thus be taken into account. 

However, equation (2) predicts that, other things being constant, 
K should be proportional to the volume of the cell. It would seem, that 
for a half cell, the value of K should then be half that of the whole 
cell. This is not borne out here, but is for the case of constriction to 
be discussed subsequently. 

The value of y/yj would be of the order of unity if we take y © 
10-* and 7 © 10-1 as above. Comparing this with the values obtained, 
we find a factor of several thousand. It might be pointed out in this 
connection that this is just the factor obtained by E. N. Harvey and 
H. Shapiro (1941) between the “relaxation” times of cells deformed 
by squeezing and cells deformed by centrifuging. The value of y for 
fertilized eggs without membranes could probably be considerably 
smaller than that estimated from unfertilized eggs. The value of the 
viscosity of dividing eggs is considerably higher than that of unfer- 
tilized eggs. Another possible explanation lies in the simplification 
of the expression for y/7 by assuming that b = 0, or that 2D; < < D.. 
It is readily shown that the effect of increasing b is to increase both 
K and y/n roughly by the factor 1/(1 — b). But this alone is hardly 
sufficient. Another factor was mentioned in a preceding sub-section 
which allows an increase of perhaps ten. These factors taken togeth- 
er could account for the discrepancy. 

An important factor not considered is that some of the para- 
meters, especially q, probably vary with time due to the numerous 
changes which are taking place. The introduction of such an effect 
should also make it possible to predict the elongation curve beyond 
the point of maximum into the region in which shortening again 
takes place. It has been suggested by R. Williamson that a factor 
which should be taken into account in this connection is the disturb- 
ance of the diffusion fields of each cell by its neighbor. This could 
result in equilibrium shapes other than a sphere. 

One might expect that »/n would be larger for cells with mem- 
branes. This is found to be the case, the increase being by a factor 
of 4.3. However, K also was found to be larger, though not by quite 
as large a factor. This might be interpreted as implying that D, is 
decreased, thus K increased, by the presence of the membrane, as the 
latter is essentially external to the cell. 

Summary. We have seen that the shape of the experimental 
elongation curve is predicted fairly well by the theory, ag illustrated 
by Figures 1 and 2. In the composite plot of values, as shown in Fig- 
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ure ii the trend may deviate somewhat at the very beginning, produc- 
ing a slight inflection. Aside from this, the entire positive part of the 
experimental curve of elongation rate with time is a parabola. The 
parameters of the elongation curves are given for various experi- 
mental conditions. The value of the parameter y/7 appears to be too 
small compared with the ratio of y to 7, estimated directly. Factors 
are discussed which can account for this apparent discrepancy. The 
parameter K is of an order of magnitude consistent with the expres- 
Sions given by equation (2). 


A COMPARISON OF THE CONSTRICTION EQUATION 
WITH EXPERIMENTAL DATA 


An expression for the constriction rate for advanced stages was 
derived by G. Young (1939b). A more detailed treatment, though 
only approximate, was given by N. Rashevsky (1941) which indicated 
that the results first mentioned might hold for a large part of the con- 
striction process. The expression is given by (Rashevsky, 1940, chap. 
iii) 


eS ee EN io (23) 
where 27 is the width of the “neck”, 7’ is an integration constant and 
Q = Kr, /88 , (24) 


in which K is given by equation (2). 

The extent to which the equation agrees with the data is illus- 
trated in Figures 2a, 2b, and 2c. In some cases the curve remained 
essentially straight or showed an inflection in about the middle part 
of the curve. Where a curve was concave downward to beyond half 
constriction, that part of the curve was used if possible. A small num- 
ber of curves could not be used, most of which were from the same 
batch of eggs. A few others were not used because the points were 
too scattered or too few. 

For each of the curves used, the point of inflection was estimated. 
Three-fourths of the curves showed inflection points below r = 26 X 
10 em. One-half of the curved showed inflection points below r = 20 
X 10+ cm., and one-fourth below 13 X 10* cm. As one measured 
values of 2r below the inflection point it became increasingly uncer- 
tain as to whether the furrow width was being measured, especially 
in eggs in which the daughter cells tended to adhere. 

The curves are not given beyond 7 since the equation necessarily 
breaks down at some value r < 7. For small widths, less than about 
6 X 10-4 cem., measurements could not be made with accuracy from the 
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photographs. Furthermore, at about this stage one would expect sec- 
ondary necking to begin. But the situation is complicated by the pres- 
ence of strands from the spindle fibers connecting the daughter cells, 

Eggs with membranes. Six eggs with membranes, two having 
been treated with KCN, were measured and the parameters 7’ and Q 
were obtained. The two KCN showed no significant difference from 
the others and were thus included to estimate the mean error. The 
average values obtained were: 7’ =1.3 + 1X 10° cm. Q=2.6 +3 X 
10-* cm.? sec.1, so that from equation (24), K =6.2 + 7 X 10% cm. 
sec.?. This value of K is somewhat less than one-half of the value 
obtained from the corresponding elongation curves. 

Eggs with membranes removed by shaking. The same group of 
eggs with membranes removed by shaking were measured for the de- 
termination of 7’ and Q, giving the following average values: 7” = 1.3 
+1X 102 em. Q=2.9 + 3 X 10-* cm. sec.+, from which K = 7.3 = 
8 X 10+ cm. sec.-?. The value obtained for this group of eggs from the 
elongation curves is about forty per cent smaller than this value of K. 
The values of Q, 7 and K are not significantly different in the eggs 
with and without membranes. An illustration of the constriction 
curve from an egg of this group is shown in Figure 2a where the 
parameters are 7’ — 1.38 X 10° cm. and Q = 3.36 X 10° cm.? sec.7. 

Eggs with membranes removed by KCl treatment. The two eggs 
with membranes removed by KCl gave the following values of the 
parameters: 7’ — 1.4 X 10° cm., Q = 3.7 X 10° cm.? sec.7, K = 8.9 X 
LO4icem. «sec. ;vand 9e Selb 10cm ae) = 4.2 <1 0 em Seca 
K = 10.6 X 10 em. sec.-*. Both these values of Q, and also K, are 
above almost all of the other values while the value of K from the 
elongation curves decreased. The significance of Q will be discussed 
subsequently. An illustrative curve from the second of these eggs is 
shown in Figure 2b. 

Eggs with membranes removed by shaking and treated with 
KCN. From these eggs the average values of the parameters are: 
rf =14+1%X 10" om. Q=—2.8 + 8 X 107 cm. sec, K = 7.2 +'8 X 
10-* cm. sec... Comparing these values with those obtained from the 
eggs similarly treated but without KCN we find them to be the same, 
that is, the curves are essentially unchanged by the KCN treatment, 
although the time to division was increased by about twenty-five per 
cent. There was no significant change in K from the elongation curves 
due to KCN treatment. 

Eggs with membranes removed by shaking, second cleavage. 
From the measurements of the constriction of the second cleavage 
eggs, the following parameters were obtained: 7’ = 1.2 + 1 X10-2 cm., 
Q=1.5 + 2 X 10 cem2 sec, K = 4.8 = 6X 10~ em. sec. Compar- 
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ing these values with those from the first cleavage we see that Q and 
K are significantly smaller here. The value of K from elongation did 
not show this change. From these comparisons we find that the con- 
striction rates are very much the same except for the half eggs. This 
suggests that the size of the cell may piay a role. But this is already 
predicted by the definition of Q and K. In Figure 2c is a sample 
curve from this group for which the parameters are 7’ = 1.24 X 102 
cm. and Q = 1.75 X 10-8 em. sec.-. 

The relationship between Q and r,. From equations (24) and 
(2) we expect that Q should vary with r,‘. But from recent studies 
(Landahl, 1942), we would expect K to be proportional to 7,2 instead 
of to 7,°. Thus Q should be proportional to 7,3. For this reason the 
values of Q from each of forty-six eggs were plotted with the initial 
radius 7,. The resulting graph is shown in Figure 3. The same sym- 
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FIGURE 3 
bols are used here as before, but in addition the values from half cells 
are designated by a plus sign. From the figure, one might suspect 
that the KCI eggs are different. In any case the correlation between 
Q and 7, , not including the half eggs or the KCl treated eggs is found 
to be 0.6 + 1. This correlation is decreased by including the eggs 
treated with KCl. The curve of Figure 3, which has but one para- 
meter, is given by Q = 0.72r.°. Thus, this prediction is at least in part 
borne out. If one uses the fourth power relation given by (2) and 
(24), Q = 220r.', the agreement is poorer. The value of Q from sec- 
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ond cleavage was close to one-half that at first cleavage. The values 
of 27, for the half cells were difficult to estimate. They were about five 
per cent less than that which would be obtained by assuming that 
each daughter cell contained one-half the original volume. In the case 
of eggs treated with KCl, the measurements could be made more ac- 
curately. The volume of the daughter cells indicated a net loss of 
about one or two per cent in volume. 

The parameter Q is roughly determined by the average slope of 
the constriction curve in the region of half constriction. If one uses 
this average to measure the constriction rate, one can say that there 
is probably a positive correlation between the rate of constriction 
and size. 

Summary. It his been found that the constriction curves are pre- 
dicted to a fair degree of approximation. From the analysis of the 
constriction curves it has been possible to predict the value of K inde- 
pendently of that from the elongation curves. In both cases the order 
of magnitude is the same, but the one set of values vary from about 
one-third to three times that of the other, without apparent correla- 
tion. Since the values of K from the constriction curves seem to vary 
with 7, as predicted, whereas this is not so in the case of the elongation 
data, one would be inclined to accept the former values in any com- 
parisons. If so, except for a possible slight effect of the KCl treat- 
ment, the only effect on K that is found is that of the size of the cell. 
The value of Q for the daughter cells is close to one-half of the value 
of Q from the original cell. 

The author is indebted to Dr. N. Rashevsky for his many sugges- 
tions, and to Dr. R. Buchsbaum for his kindness in permitting the use 
of the photographs and records. 
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TABLE OF PARAMETERS 
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Elongation Curves 
A sec-t .025 024 .020 .020 .025 
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Constriction Curves 
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_ When total cell number is used as the basic parameter of growth, 
rational equations which describe colonial and organism growth under 
varying circumstances have been derived from a single differential form. 
These equations result from making specific, but reasonable assumptions 
about two additive factors, ¢ and 6 which determine community growth. 
The first factor (¢) is assumed to arise from conditions within the grow- 
ing cell itself, while the second factor (6) arises from interactions be- 
tween the growing cells of the community. If it is further assumed that 
the cells of a community are homogeneous with respect to density and 
volume, it has been shown that the mathematical expressions commonly 
used to describe growth data may be rationally derived from the general 


form. 


In attempting to describe rationally how a cell colony or a multi- 
cellular organism changes in time (t), the crucial preliminary ques- 
tion is to decide upon a parameter of description, y(t). The necessary 
attributes of this parameter are: (1) w must, as nearly as possible, 
depend only upon the activity of living cells. (2) w must be such that 
all other parameters can be unambiguously expressed in terms of y. 
In view of these requirements it is strange that most theoretical stud- 
ies of community growth have chosen for y virtually every dependent 
variable except total cell number. (See, for example, Ludwig, 1929; 
Bertalanffy, 1934, 1938; Wetzel, 1937). It will be shown that the total 
cell number, N , can best satisfy the conditions on a growth parameter, 
and furthermore, that in terms of N all the important rational growth 
equations and two of the equations heretofore regarded as empirical 
can be shown to be but special cases of a fundamental differential 
form. 

The phrase used above, “activity of living cells,” is a vague term 
which veils much ignorance both about “living” and about “cells.” It 
is not the place here to investigate the term. Let us only remark that 
Fen eee ee eiag, Ontidal Projeel NG. BEL ORE, 
Unit A-8. 
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within the limits of our knowledge a living system always exhibits 
the phenomena of metabolism and irritability. Under this definition 
not all the substance of a cell is living, since yolk stores, fat stores, or 
water stores, within a cell are no more “living” than the same sub- 
stances outside the cell. Nevertheless as the absolute amounts of these 
substances in a cell change with time, the usual parameters, total 
mass, surface, and volume may also change. It is thus very clear that 
the parameters, mass, surface, and volume could not possible be ex- 
pressions—except in a very indirect sense—of the phenomena we call 
living. Furthermore, indices of metabolism and irritability vary in 
time as a result of physiological processes not associated with growth, 
e.g., exercise, temperature changes, etc. On the other hand the pro- 
cess of cell division is unquestionably associated with living proto- 
plasm in a most direct way, so that it must be concluded that the 
measure of cell division, namely cell number, is a basic parameter for 
a theoretical analysis of growth. 

To derive the fundamental differential equation describing 
growth, consider a community, be it a colony or an organism, of N 
cells at time ¢. Let us postulate that two factors are responsible for 
changes in N as t increases: (1) a factor ¢, which resides within the 
cells, and would be present were the cells infinitely removed from 
each other, and (2) a factor 6, affecting the cells because of interac- 
tions between cells. Assuming that these two effects are additive, we 
may write, 

dN 

ra i Poe wlan (1) 
(These are not stringent postulates; we only require that both cell and 
interaction contribute to the division of the cell, as of course they do.) 
Let us now consider ¢ and @ in turn. Assuming no interactions be- 
tween cells, it is reasonable to suppose that at any one time the rate 
of cell number increase is proportional to the number of cells pres- 
ent. However, as time progresses the “constant” of proportionality 
will in general change; in other words the “constant” of propor- 
tionality should be a function of time. Such an assumption would lead 
to the general form of differential equation, 


dN 
FE 1 (rs Besos Un, t) 

where the «’s are physical parameters and t is time. This is indeed 
what we are about to write. While this equation is mathematically 
quite permissible, it obscures an important step in physical reason- 
ing, namely, the way in which ¢ is introduced into the right-hand 
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member. For the pointing out of this step, and for the following in- 
clusion, we are grateful to N. Rashevsky. Physical functions are de- 
termined at any instant by the values of physical variables, and de- 
pend on time only insofar as the latter depend on time. Thus it may 
be that a parameter, x; , changes in time, 


dx; 
FE Ie Tans Pm) 


whence, 
xi = fs (%i,, t) 


where «;, is the value of x; at the zero point of time. Substituting in 
the equation for N, we are then able to replace x; by x:,, t, and simi- 
larly for other parameters. The right hand member then becomes rig- 
orously meaningful, with the zero point of time operationally defined, 
and can be evaluated for any given t. 

Furthermore we shall later cite cases where the coefficient of N 
may depend upon N as well as t. Accordingly we have in general, 


d6=A(N,t)-N. (2) 


On the other hand, while the interaction effect must somehow depend 
upon the number of cells present, the dependence need not be linear, 
but rather upon N*, where in the simplest case k is some constant. As 
before, the coefficient of proportionality is a function of N as well as 
of t. Thus. 


6=B(N,t)-N*. (3) 
Putting together (1) (2) and (3) we have 
NAAN tN BAN sh) GND (4) 


dt 


The fundamental character of equation (4) is demonstrated by the 
fact that by making reasonable assumptions about the functions A 
and B, many phenomena which affect growth can be incorporated 
into the form (4) without sacrificing its rationale. 

The intrinsic factor ¢. For instance, the phenomenon of the lag 
phase in bacterial and protozoan growth can reasonably be regarded 
as the approaching of the steady state initial to active proliferation. 
In the simplest case, the effect might be expressed by the well known 
transient function, & e*’, in which case we would set 


Ae ee (CUA). 


Clearly as t > », A > C, = constant, but for t very small, i.e., during 
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the early stages of growth, A might be considerably below its eventual 
value, and 4 correspondingly low. 

The interaction factor 6. It is evident that 6 will express the net 
effect of inhibitory and facilitative interactions. If these interaction 
effects are independent of one another, each will contribute an addi- 
tive term to 6, and 9 will be represented by a sum; if they are all re- 
lated, 6 could be a single product function. In the most general case, 0 
would be a sum of product functions, 6; , where each 6; is given by 


6,=Bi(N ,t) -N*, (5) 
and 
=> 6;. 


As examples of negative or inhibitory interactions (6; < 0) we might 
cite (1) the competition by cells for nutrient materials; (2) mechani- 
cal effects, e.g., compression; (3) spatial effects, such as the distor- 
tion of nutritive diffusion fields of a cell by its neighbors; (4) produc- 
tion of harmful metabolites. Consideration of positive interactions 
immediately brings to mind endocrine effects, e.g., “growth” hormone, 
thyroxine, etc. One way of treating these effects would be to consider 
a gland as a sub-community which does not exert its influence on the 
entire community until some time characteristic of the gland. The 
simplest mathematical formulation would then be 


6; = B; (function of + only) - N;?, 


where B must express the growth curve of the gland alone, and N; the 
number of secreting cells in the gland (which are implicitly assumed 
functionally homogeneous). An alternative and simpler way is, how- 
ever, suggested by the experimental fact that the most important case 
of hormonal] influence, as, for example, puberty in mammals, occurs 
at the same time that d?N/dt? = 0. It would thus be equally reason- 
able to set 

d?N 

dt? 


Buf 


and therefore 


@N 
== Se, eee 
8:= Nef. 


This is an important case where B is a function of N as well as of t. 
The exponent k; (5) may be chosen to conform with any particular 
hypothesis. Concrete illustrations will be presented below. 

These examples suffice to show that under certain specific hypoth- 
eses we can express interactions quite rationally by the function (5), 
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and thus complete our justification for the use of the form (1). We 
will now study integrals of the form (1) under various rational as- 
sumptions regarding ¢ and 6. 

Case I: Linear Proliferation; Universally Deleterious Metabolite 
Produced by all Cells at a Uniform Rate. If we assume that cells in- 
finitely removed from each other have equal and constant division 
times, it is easy to show that, 


@—AN, (6) 


where A is a true constant (no transient or lag effects). If now the 
cells are brought together and if each cell produces a substance which 
is harmful to every cell in the community including itself, and the 
proportionality constant is B , then the sum of harmful effects is BN?. 
This follows because the total number of effects of N cells upon them- 
selves is BN; Furthermore each cell acts upon N—1 others, and there- 
fore the number of effects of cells upon their fellows is BN(N-1). 
The total number of harmful effects is thus BN + BN(N-1) = BN?, 
and we may set 


6= —BN?. (7) 
Combining (1), (6), and (7), we obtain 
dN A 
—= — 2— BN (= — : 8 
dt AN — BN?=—8 v7 N) (8) 


If B =f and A/B =a, (8) integrates directly to 


a 
Sa ae ee a Or ee 9 
N= T+ Ke “ 


where K is a constant of integration. 


Case II: Linear Proliferation; Universally Deleterious Metabolite 
Produced by all Cells at a Variable Rate. In addition to the hyptheses 
of Case I, let us postulate that the proportionality constant B changes 
in time as discussed earlier. We can attempt to describe this time 
variation by a power function of ¢, thus 


B=> Rpt. 


8=0 
Substituing into (1) we obtain, 


T= AN - MS It. (10) 


s=0 
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Equation (10) integrates to 


a 


N = const. + uieee 


(11) 
Case III: Constant Rate of Proliferation; Foodstuffs as Limiting 


Factor: An interesting simple case is to consider the proliferation 
rate, ¢, (1) a constant, 


o=ak (a, k, constants). 


(Here we have broken up the constant into two factors for algebraic 
convenience), and to assume that retardation is due to a shortage of 
foodstuffs. The latter effect can evidently be set proportional to the 
number of cells present, let us say, 


CNA 
The rate equation, 
aN 
— = ak — TZ 
ete EN, (12) 
this integrates to 
N=, Glee se mt e (13) 


Case IV: All Cells Produce A Stimulating Substance Which takes 
Effect According to Pharmacological Laws; Foodstuffs as Limiting 
Factor. If all cells produce a growth stimulant, the amount of sub- 
stance present at any time is proportional to the number of cells. The 
effect, however, need not be so; in fact it is common pharmacological 
experience that the effect of a dose is proportional to the per cent in- 
crease in the substance. (This general principle has much wider ap- 
plication than pharmacology; in physiology it has been formulated 
quantitatively by Weber and Fechner). Under those circumstances 
the effect of unit amount of substance would be proportional to the 
logarithm of the amount of substance, which in turn is proportional 
to the number of cells. The total effect is then proportional to the 
effect of unit amount and to the amount, in other words, 


¢=—AaN log k.N 
as before, 


6=-— const. xX N=— (alogk,)N. 
Combining by (1), 


aN 
GE ON log kN — (alog ky) N. (14) 
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Equation (14) integrates to 


‘A 


ky ceat 
N= a € (15) 
where c is a constant of integration. 

These four examples we believe are sufficient to indicate the use- 
ful generality of the form (1), which in turn was based on the use of 
N as the essential parameter of growth. Before discussing them fur- 
ther we will consider the relation between the differential and inte- 
grated functions (8)-(15) and their analogues using mass, M , as the 
essential parameter. Let us begin by assuming that for the k-th cell 
the average density is j, and the volume, v; , so that its mass is 


Mx = pKVk . (16) 
and obviously, 
k=N 
M=>=> pPrVx- (17) 
k=1 


Since M =f (N, t, fx, vx) we may also write, 
aM OoMdN **oM d px ken OM. dv; 
= +3——+>— 


aie iranian ae (18) 
dt ON-dt and dt ga oy dé 


The set, (17) and (18), obviously give the desired relations. However, 
it will be observed that in practice an inescapable assumption to be 
made is that the cells among each other are homogeneous with respect 
to mass, density, and volume. (Although our differential form did not 
depend upon any such assumption, each of the specific hypotheses J to 
IV contained it.) 

In that case (18) becomes simply 


M=N pv. (19) 
Whence on substitution into (18) we have 


dM MdN pace eee ae 
Gi Nd pdt vdt oo 
If in addition, p and v are considered constant in time, (20) becomes 
simply dM/dt = const. X dN/dt , and the four cases we have studied 
can be written in terms of mass and time, using the same contants 
referred to mass and time. The mass analogues of Equations (9), 
(11), (18) and (15) are thus; 
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Qa 

=o (21) 
oa ah Kea 

= —_. 22 
M=const. + aaa (22) 
Mai =e); (23) 

and 

UE gooe, (24) 


iD 


which may be recognized respectively as the equations of T. B. Robert- 
son (1926), R. Pearl (1925), Mitscherlich and others such as R. M. 
Jenss and N. Bayley (1939), and B. Gompertz (1825). The first two 
of these have heretofore been derived on the basis of chemical anal- 
ogies which were admittedly fictitious while the last two have been up 
to now regarded as purely empirical (see A. J. Lotka (1925), S. Brody 
(1927), S. A. Courtis (1932) and I. Schmalhausen (1931). 
In closing we may remark, perhaps unnecessarily, that there is 

a many-to-one correspondence between an integral form and a differ- 
ential equation; so that interpretation of the constants in equations 
(21) to (24) is not unique, although the meanings ascribed to them in 
our derivations are certainly one possible set. 

We are grateful to Professor S. F. Cook and Mr. F. Kreutzer for 
reading and criticizing the manuscript. 
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A NOTE ON THE DIFFUSION OF ELECTROLYTES IN CELLS 


ALSTON 8. HOUSEHOLDER AND ROBERT R. WILLIAMSON 
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A somewhat simpler solution is given to the problem previously dis- 
cussed by R. R. Williamson relating to the diffusion of a ‘metabolized 
electrolyte whose ions have different mobilities. 


R. R. Williamson (1939, 1941) has calculated the electrical forces 
that would be present in a cell metabolizing electrolytes, making cer- 
tain simplifying assumptions and approximations. The purpose of 
this note is to introduce some substitutions which simplify the form 
of the result, and to correct an error in the second paper. In the first 
paper the problem considered is to determine the electrical forces un- 
der the assumption of differential transport rates for the two ions 
within the cell, but infinite permeabilities and external diffusion co- 
efficients. Williamson’s fundamental equations (9) whose solution 
gives the “average” internal concentrations of the two ions within 
the cell are derived from a pair (not numbered) which may be writ- 
ten 


Ts dD. & 
Pu ¢,? = PY CrGe ae 6D.¢, gee Mest O es “ee 6D.Co. = 0, 
: ‘ 1 
ao ee (1) 
Par EF p Mo + 6D_¢_ — rq. — 6D_¢)-=0. 


The concentrations c are all given in gm cm’, and the productions q 
are in corresponding units. Evidently mol cm is a preferable unit, 
especially since the molar production rates, and the external molar 
concentrations are equal. Moreover, the interest lies in the difference 
between the internal ionic concentrations, and this difference can be 
non-zero only so long as the production rate is non-null and the two 
diffusion coefficients are different. It is therefore appropriate to set 


_ Cor Co- fe as a 9 
Oe Ym, MALE get nny ah Me (2) 
to obtain, after an obvious reduction, 
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GVO" GR. Cre eC ane 6 
(iz) . Mise OS 

(3) 
CNS? CoC One eee Coe 6 ; 

= jee ae —-—— — — g — —¢, — 0. 
(iz) MiM ius Mam es 
Since for D, = D_ we have é. : M. =é_: M_, the suggestion is immedi- 
ate that we introduce parameters which vanish with the differences. 
We have, in fact, on adding and subtracting the equations, the equiva- 


lent pair 
Cy CoN Cx is 6 Tsai 1 
Ye Se 
(az, ral -) p ios =) 


5 ‘ e ; (4) 
Cs = MA TGd i .Cz é..\ (ie “5 = 
Gin ae) +7 Get ar) = > 8 to 
If we set 
psec, Ga as a 
ii Pas) oe Lage ge 
ae Fo 1 linet iog it ws 1 412, 4 36 (5) 
ye p? Di a v p & B: “das 
we have 
rhe ep 
xy ene He 
(6) 
etl (yt) at, 
p p 
which, for small wu, in fact for 
wv — ey 
p 
can have but a single solution. We have thus 
v= 27F-u=0 (7) 


as a fairly simple cubic for determining the concentration difference 
for the two ions. Being already in “reduced” form, the closed solution 
is easily written. However, the expansion in powers of x contains 
only odd powers, and is given as follows: 


xe=u+ vu + 3v7u5 +.--, (8) 
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For a cell with finite permeability (Williamson 1941), if we set 


— mY? N &? 
esi Ta 


R. R. Williamson’s equations (12)-(15) require correction and should 
be as follows: 


Qn? = 
6D, Cr Cu t SAC., (9) 
Qro 6 
Bh Ow Cot B= (A + A) (Cur +O), (10) 
SE Fics kn BAC 
ap. C7 Gx — EAC-, (11) 
Che (At ACC 
ee 1— 0 To 1) ( i- o) > (12) 


where now the C’s and the Q are in terms of mols per cc. Then on 
adding (9) and (10), (11) and (12), and taking the sum and differ- 
ence of the resulting equations we obtain the following equations: 


05 re Toe 4 Tae Mt ite = 


Ee aoe Sho oh: 


is = Bs 3 
Ge ae C_.— pd ta EA? == es (A I A,) A, 
0 
To To” ae een vals 
Q (60, 6D. 3h, 3h. ) 


AP EAC, #C) +E5 (A + A,) (C. +. C. + :26,). 


If we introduce the following notation: 


Baral tats pact tabs 
Mie eae eee ref 


EA tae Te 
10=0| ay 6D. Bh, a | 


e (13) 
deme] 8 t,=A:, 
y=C,+C_, Y,=Cy+C., 


E, = £6/"%. 
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Our equations are 


AQ=y — 20, HER? PEL Fi), (14) 

AC=at Exy + &,(" + 41) (y + 2Cd). (15) 
Since 

eae oo 108s 
Yo To 
and 
| (% + 21) %1| © x? 

usually ; 


a2 >> b: (et OD Gare SE 0, 


and our equations are very closely approximated by 


AQ SY Ce et a (16) 
AQ=2+ Exy. (17) 

From (16) 
Y—AQ: -+- 205 6e-, (18) 


and substituting in (17), 
AQ =e sant 2Cye2er le 


or, after rearranging, 
&*a AQ . 
$$ * 7 t+ ————_— = 0 
1 +&(AQ + 2C,) 1+ €UQ + 2C,) 


If then we let © 
oe ie AQ \ 5 
1+ EQ 20a y a) 1 BE OS eC 
the equation has the form of (7) and the expansion (8) is the same 


when subscripts are placed on the wu and v. The expression for y is 
then obtained by substitution into (18). 


This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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_The periodic and resonant properties of a closed neuron circuit are 
exhibited and applications to some visual and electrophysiological phe- 
nomena are discussed. 


Various properties of closed neuron circuits have been discussed 
previously. In particular, Landahl and Householder (1939) discussed 
the conditions for periodicities in the self-exciting neuron. 

In the circuit to be considered, neuron J, producing excitatory 
substance « only, synapses with JI, which produces inhibitory sub- 
stance 7 only, and JJ in turn synapses with J. External stimulation of 
amounts #, and H, may be introduced through fibers I’ and II’ at 
synapses (1, 2) and (2,1) respectively. 

We consider the equations 


d é 

Fa bs (Ba — §— In) — He , 

mn (1) 
a Ps (Bs each easy 


where a; and f; are reaction constants having positive values and E; 
are arbitrary functions of time. The auxiliary equation for this sys- 
tem has the roots 


EIAs Cas 2 Gp HaOas\ BF ah 9 
a eo: a e ‘ (Par ) 01102 — fi Bo (2) 

The behavior of the system will be periodic provided 
(A) — ee — Be < 0. (3) 


It will henceforth be assumed that this condition is satisfied. 
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We will adopt the notations 


0, + Og == 27, 
01, Og + Bip2=g ’ (4) 
p—g=-@. 


The complementary functions may then be written 
&- = EF (c,e'4' + e744) , 


= a ; , ; t (5) 
95 —€ (6,0 One 2), 
or, by using trigonometric functions, as 
&—e"'(A cosqt+ Bsin qe) , 
(6) 


je— e*' (A’ cos gt + B’ sin qt) . 


The constants A’ and B’ are not independent of A and B, but satisfy 
the relations, 


Oz — a (7) 
B=|aA sy lop on 
Making these substitutions and rewriting, we have 
E> — e*' » sin (qt + ¢) 
a (8) 
7, = e% Vz y sin (gi + +0) 
1 
where y and ¢ are the constants of integration and 6 is given by 
en ice ee og ee (9) 
V BBs 2V Ps 


In the case that the imposed stimulations, HE, and E,, are con- 
stant, the particular integrals for « and j are given by 


Em = Bilas (H, — hy) — B.(H, — he) 1/9, 


: 10 
Jn = Bala (Ea — he) + B,(Ea — ba) J/g a 
The complete solutions for constant H, and E, are then 
e=e"»sin(qt+¢) + & 
(11) 


fle 
ee Fr sincat +6+6) + Jo. 
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Since « and 7 must be always positive we must satisfy the conditions 
2 0. (eens (12) 


This is equivalent to saying that in the E,, Ez plane, the point (£,, E2) 
lies above both of a certain pair of lines 


de (E, — hi) — B.(E, — ho) =0, 
Bi (BH, — h,) + 0, (F, — he) =0. 
We see that we must have 
Le ies 


(13) 


(14) 


In evaluating y and $ we will assume that at t, = 0, 
Big alts 2d 5 5 Ol; E,-h,+e 20, (15) 


with #, and E, satisfying conditions (14). Evaluating ¢ using these 
restrictions, we obtain 


&=ysind + &, (16) 
= =p, (he — Ra 4y). ato = Hy COS oO — Py Sine. LY) 
dt t=0 

Then 
q (Eo =I #5) 
SSS (18) 
CD sei He 2) aS Da Fe) 


+ 
y =A lA EB: Seyi Oren Dene Oe le O° (eo ca)*| . (19) 
qd 
In order for the function e(t) to be continuous for all positive 
values of t we must have « = 0 when de/dt = 0 and d?</dt? > 0, or 
&, + ery sin(qt, + 4) >0, (20) 


where t, is the time from t = 0 to the first minimum of ¢. We evaluate 
t, by setting the time derivative of ¢ equal to zero and find 


t,=* [ sinx(—2) - : |. (21) 

qd Vg J 
Corresponding conditions are obtained from the requirement that 
7 2 0 at its first minimum. 
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The case of sinusoidal applied stimulation. 
We will consider the case in which 


E,(t) =E, sin 6¢ + Ee, 


(22) 
Ei, (t) =yE, sin 6t + whe. 
For a particular integral of (1) we choose 
V=c, sin (6¢ + w) +c cos (66 +m) +63. (23) 
Equating coefficients, we find 
= (@? es 9) KE» 
“(2 = g)? + 4p?’ 
2p0KE, 
>= =, 24 
(2 — 9)? + 4p6? aa 
Ca Bil (de = fou) BE» = Oeltg ++ Bshez] /9 = Ew » 
where 
Ke == 8? [| (Ge =few) = Bees, (25) 
and » and g are given by (4) above. 
The particular integral for « becomes 
spas = sin tee @ ed) en (26) 
p 
where 
= 1 6? ‘ eee on 
ie mie g2 + 2 a+ 0-2)" |, (27) 
with 
—_— Qe — Bou F (28) 


The maxima and minima of p(@) are found by setting dp/dé = 0 
and solving for 6. This gives 


Perens ase) ee dg) ee 0 (29) 


for the extreme values of p(@). 
The amplitude at 6= 0 is given by 


Ey 
i = p,rE,/g . (30) 


For very large 6, 1/p, and hence the amplitude, approaches zero. 
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Discussion. 

To recapitulate, the above formulation shows that when a closed 
circuit consisting of one excitatory and one inhibitory neuron is sub- 
jected to appropriate constant stimulation an oscillatory response is 
set up having a period characteristic of the circuit, the oscillations 
diminishing toward a steady value characteristic of the stimulus. It 
is also shown that when a sinusoidally varied stimulus is applied the 
circuit responds synchronously, the amplitude of response being de- 
termined by the frequency of the imposed stimulation. The frequency 
to which the circuit responds maximally is known as the resonance 
frequency, and is related to the natural frequency of response to 
steady stimulation. 

There are many instances of animal behavior under natural and 
experimental conditions which have the properties described. In some 
cases there is evidence that the above described mechanism cannot be 
considered as operative, at least exclusively. Of these we may men- 
tion the rhythmic potentials recorded from the brain, which have been 
shown by Libet and Gerard (1938) and others to depend on other 
factors than anatomical relationships of neurons. Such mechanisms as 
the scratch reflex and various cases of oscillatory (clonic) response to 
a steady stimulus would conform to the scheme presented above, with 
the proprioceptors, especially stretch receptors, corresponding to the 
inhibitory element and the motor nerve corresponding to the excita- 
tory element. The above formulation would be adequate only within 
such an intensity range that the linear assumptions could be con- 
sidered to hold. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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A THEORY OF ELECTRICAL POLARITY IN CELLS 


ROBERT R. WILLIAMSON AND INGRAM BLOCH 
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In general, there exist electrical potential differences along the axis 
of a spherical cell containing an asymmetric diffusion field. The elec- 
trical potential is calculated for any point on the axis of such a cell, 
ae is assumed to contain metabolite ions and charged colloidal par- 
ticles. 


Consider a spherical cell of radius 7,, with its center at the 
origin. We shall assume that an existing metabolic polarity (Rashev- 
sky, 1940) makes the concentration gradients greatest along the 
x-axis, and treat the cell by halves: the first half is the hemisphere 
whose axis is the positive z-axis, while the second half is the hemi- 
sphere whose axis is the negative z-axis. Our problem will be to de- 
termine, approximately, the electrical potential at a point P on the 
x-axis where x=ar,, if —1 =a = +1. Thus we shall be able to 
determine the potential difference between any two points on the 
x-axis. 

We shall deal with + and — ions which are being produced or 
consumed in the cell, and with catalyst colloidal particles which are 
assumed to be present in the cell. For simplicity we assume the perme- 
ability of the cell wall to the ions is «, and the external diffusion 
coefficients are also oo. 


Cj. = average concentration of ion in moles cm- in the 7-th 
half cell 

Q; =average rate of production of ion in moles cm-* sec” in 
the colloidal charge in the i-th half cell 

D. = internal diffusion coefficient of ion in cm? sec 

= average difference in concentration of ions of opposite 

signs in moles cm- in i-th half-cell 

B. =wmobility of ion = D./kT 

k =Boltzmann constant 

T =absolute temperature 

Is 

fis 


BS 


4 


= transport of ion in moles cm” sec* 
= electric force on a + ion in 7-th half cell (fi: = —fi-) 
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—, f i+ To 

DE, 

V’; = electrostatic potential in the first half cell resulting from 
the metabolite which is in the 7-th half cell 

V"; = electrostatic potential in the second half cell resulting 

from the metabolite which is in the 7-th half cell 


n; == number density of charged colloidal particles, inhibiting 
Q in the i-th half cell 
fila TP 105 
eee 
y =number of excess + charges in units of one electronic 


charge, on a catalyst particle 
é =charge on an electron = 4.8 X 10° E.S.U. 
6; =%nve=charge density in i-th half resulting from charge 
on colloid 
é; — electrostatic potential in the first half cell resulting from 
the colloidal charge in the i-th half cell 
¢’; = electrostatic potential in the second half cell resulting 
from the colloidal charge in the ?7-th half cell 
pi =V', + ¢:=total electrostatic potential in the first half 
cell due to the charge in the 7-th half cell 
yp", =V"; +6"; = total electrostatic potential in the second half 
cell due to the charge in the i-th half cell 
N = Avogadro’s number 
RK =+Vy? + 2 
p =charge density 
Our steady-state diffusion equations become (Rashevsky, 1940; 
Williamson, 1939, 1941) 


270Q: 
3D, = (54-2) C26. 4c; 
277Qi 
3D. = (5 oa F,) C,- te C.- a ACs 
27.°Qs (1) 
3D, aoa (5 a F,) Co. aa Cr a AC 
21?Q2 
SD (5 — Fi) Ge C4, 


If we can determine F’, and F, , we can solve these four equations 
for the four unknowns, C,., C:., Co, Co. A knowledge of these, plus 
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a knowledge of the quantities 6; resulting from the distribution of 
charged colloidal particles, will enable us to calculate the approximate 
electrical potential as a point function, in particular, at the point 
FAS =a, y—0,2=0) 


ees 
mer: ; Toy 
where 7 is the volume of the cell and K is the dielectric constant 
p=NeA,+6, ifa>0 
SSN GARE het 
dr = 2n RdRdx 
r = distance from dz to P 


72 — Re? SE (ar, = wee 
Then 


, _ 2n(NedA, + 6) (er a (pie RdR 
alee K 0 VR? + (a% — x)? 
Assuming a > 0, we find 


-  Qnr? (Ned, + 6:) 3a 
: Sy ae 1 (a +1) | , 


Me 3Ka 

Qnro?(NeA, + 32) f , , 3a 
————————— —— ef (1 2) 8/2 
y's 30K | « Se 5 Tea ea) 


y (aro) = W's oe Y's 


__ 2ar,? ia : 30, : ae 
= Bika | Ne e ea ere ( a) | 
3 (2) 
e a 
—(Ned, + 4,)| 2a" -~S+1- (4 a)” || . 
Soe Oe ey () 
be 0 EY dar) ST do’ 
2ary2vE “ (1 — 20?) (1 +. a2)? —1 
Fee 2 2) | 2a - AA 
yi Tae Oe On) a a | 
(4) 


ee le 


a? 


— (Ned, + 6,) [40 5% 
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Since in performing our integration we assumed « > 0 the foregoing 
computations apply only to the first half of the cell. The results may 
be applied to the other half of the cell, however, if we exchange 
NeA, + 6, and NeA, + 6. Hence for any two points on the z-axis 
equidistant from the origin (one in each half) defined by x = = a7 


y (ato) — yp" (aro) 


Qathig* bs . 
= re |tedteas +8:) — (Ned, + ) [2a Pye oe 2. (5) 
21 97 EV e (1 — 2a?) (1 + a2)? -—1 


ea acie 5) > 


Qe 


= NeAseeas) [40 in 


F,and F, , expressed as the difference of two small quantities, in them- 
selves small differences, will be much less than 5 in absolute value. 
Hence, if we let 


and neglect F in (1) we have 


Cy. + Cy. — Co, — Co- = C, — C, = 6(Q; — Q.) () 
and 
aod ne eto. A 
A, = Cy. - Cy = f+, 
DLE aie 
5 5 
or 
A, — 42 = €(Q, — Q2). (8) 


Now let us assume that all the colloidal particles are positively 
charged. Then, treating them as a gas in a field of force derived from 
a potential, the distribution of the particles will be 


nN Ey (UONeLe » ( To ” 
cabernet ea (£2) -y (3) ] (9 
Ny ) 
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where y = 8Nv/2, v being the volume of a colloidal particle. 


Let 
Y= 1, — Ne 
Q,=Q) — an, 
Q2=Qo — an, ’ 


where a is a constant. Subtracting the last pair 
Q: — Q2= —a(n, — ne). 


If we remember that 6, = nye we can obtain the exponent of (9) asa 
constant times y 


‘Spit aeS OnE = ak vy (To/2) — yp" (— re/2) | 


= Y cat re [Netw as ve: (10/5 — 19) | te) 


6KkT 
or 
Ne —— ely 
Ne 
where 
6 VE 
u=frri(e-a) 
(10) 
ze ver a 
p= (10\/5 — 19) — 101/5 = 197 
By the same argument as that in N. Rashevsky, 1940, we find 
M 
y =2n tanh =. (11) 


A stable distribution of colloidal particles is possible if equation (14) 
has a real positive root. If that root exists and if My is very small, 
we may expand tanh My/2 and keep only the lowest powers. 

We get 


M?y? 
Be TM (1 PTE 
whose roots ~ 0 are 


ne to a 8 (12) 


assuming nM > 1. 
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We have now obtained all the equations necessary to a solution 
of the problem, and we can predict the magnitude of the polarity and 
electrostatic potentials. In a subsequent paper, we will discuss the 
conditions for polarity, the kinds of polarity that can result from this 
mechanism, changes of polar state of a cell, and orders of magnitude 
of the electric potential differences at the poles. 

The authors wish to acknowledge their indebtedness to Dr. N. 
Rashevsky and Dr. A. S. Householder for their frequent and helpful 
criticisms and suggestions. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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In connection with a previous paper on the same subject this differ- 
ent approach to the problem is outlined, which also leads to an expression 
for Yaa values of different binary combinations of musical 
sounds. 


In a previous paper (Rashevsky, 1942) we have treated the pres- 
ent problem under the assumption that the intensity of excitation of 
the centers corresponding to overlapping harmonics was the same as 
for centers corresponding to non-overlapping ones. We shall treat 
here the problem under a somewhat more general assumption. We 
shall consider the intensities of excitation of the centers correspond- 
ing to the non-overlapping harmonics as the same, while for a center 
corresponding to two overlapping harmonics, we shall take the exci- 
tation to be twice as great. 

The total intensity of excitation of n’ mutually inhibiting centers, 
excited with different intensities e; , is (Rashevsky, 1942a) using the 
same notations as before 


eae (Ne 1) | ihe. . (1) 
where n; is the number of centers with the intensity of excitation e; . 
We have in general 
Sng 1. (2) 
In our particular case we have 
Ga 6 ed Na ae 16g he (3) 
Hence equation (1) becomes 
E=[1—-— (v —1)6](2n, +7). (4) 


Introducing as in loc. cit. the quantity x , we have: 
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van; n= /e=—n'p/ (2p —-1); 
n',=n/p=n/ (2p — 1); (5) 
n',=n — n'/(2p — 1) = (2p — 2)/(2p — 1). 


Hence 

Qn', + n',=2n'/z. (6) 
Therefore, neglecting 1 as compared with n’, and putting n = an, 
bn = B, we have 

E=2n(1 Baz): (7) 


E decreases monotonically with «. Since experimentally it is known 
that the pleasantness rating has a maximum for approximately « = 1-8 
(Rashevsky, 1942), therefore we cannot in this case identify E di- 
rectly with that pleasantness rating. We may consider however, that 
the pleasantness rating is determined by the excitation intensity of a 
higher center, which is connected with the center at which the exci- 
tation is H , by means of neurons having the characteristics: 


Ag yc 

bs 
and for which E serves as a stimulus. For such a neuron in a steady 
state the intensity of excitation E* will have a maximum for a certain 
value of E' (Rashevsky, 1938). By a proper choice of the parameters 
involved in (7) and (8) we may obtain approximately the desired 
rank order of pleasantness. The difficulty with the octave, found in 
loc. cit. vanishes in this treatment. The octave will have a higher rat- 


ing than the fifth, because for the octave x = 1-5, as against 1-67 for 
the fifth. 


ASS BES a Sb: (8) 
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CONSERVATION OF SCHOLARLY JOURNALS 


AN ANNOUNCEMENT 


The American Library Association created this last year the 
Committee on Aid to Libraries in War Areas, headed by John R. Rus- 
sell, the Librarian of the University of Rochester. The Committee is 
faced with numerous serious problems and hopes that American schol- 
ars and scientists will be of considerable aid in the solution of one of 
these problems. 

One of the most difficult tasks in library reconstruction after the 
first World War was that of completing foreign institutional sets of 
American scholarly, scientific, and technical periodicals. The attempt 
to avoid a duplication of that situation is now the concern of the Com- 
mittee. 

Many sets of journals will be broken by the financial inability of 
the institutions to renew subscriptions. As far as possible they will 
be completed from a stock of periodicals being purchased by the Com- 
mittee. Many more will have been broken through mail difficulties 
and loss of shipments, while still other sets will have disappeared in 
the destruction of libraries. The size of the eventual demand is im- 
possible to estimate, but requests received by the Committee already 
give evidence that it will be enormous. 

With an imminent paper shortage attempts are being made to 
collect old periodicals for pulp. Fearing this possible reduction in the 
already limited supply of scholarly and scientific journals, the Com- 
mittee hopes to enlist the cooperation of subscribers to this journal in 
preventing the sacrifice of this type of material to the pulp demand. 
It is scarcely necessary to mention the appreciation of foreign insti- 
tutions and scholars for this activity. 

Questions concerning the project or concerning the value of par- 
ticular periodicals to the project should be directed to Wayne M. Hart- 
well, Executive Assistant to the Committee on Aid to Libraries in 
War Areas, Rush Rhees Library, University of Rochester, Rochester, 
New York. 
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SCOPE OF THE BULLETIN 

1. The Bulletin is devoted to publications of research in Mathe- 
matical Biophysics, as contributing to the physicomathematical foun- 
dations of biology in their most general scope. 

2. Papers published in the Bulletin cover physicomathematical 
theories as well as any other mathematical treatments of biological 
phenomena, with the exception of purely statistical studies. 

3. Mathematical studies in physics or in borderline fields in 
which a direct connection with biological problems is pointed out are 
also accepted. 

4, Emphasis is put upon the mathematical developments, but a 
description and discussion of experimental work falls also within the 
scope of the Bulletin provided that description or discussion is made 
in close connection with mathematical developments contained in the 
same paper. 

5. Outside of the scope of the journal are papers of purely sta- 
tistical nature or papers concerned only with empirical equations. 


PREPARATION OF MANUSCRIPTS 

All manuscripts should be typewritten double space. Equations 
should be numbered on the right consecutively. Do not use such nota- 
tions as “equation 2a” or, “equation 5’” etc. References should be all 
given at the end of the paper, arranged in alphabetic order by authors 
and for each author chronologically, following exactly the style used in 
the Bulletin. In text, reference should be made by giving in paren- 
theses the name of the author followed by the year of publication. In 
case of several publications by the same author in the same year, use 
notations “1940a’, ‘1940b’, etc. 

In writing equations, slanted lines should be used wherever pos- 
sible. 

Every paper is to be preceded by a short abstract. 

Drawings should be prepared in a professional manner on white 
paper or tracing cloth following as closely as possible the style used 
in the Bulletin. They should be drawn approximately twice the scale 
of the finished reproduction. Lettering should be made by means of a 
lettering guide. If desired by the authors, drawings may be made at 
the Editorial Office according to author’s instructions, at cost. 
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